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Abstract 

The decays of XcJ tt^tt^ , (J = 0,2) are discussed within the standard and modified 

hard scattering approach when including the contributions from twist-3 distribution amplitudes 
and wave functions of the light pseudoscalar meson. A model for twist-2 and twist-3 distribution 
amplitudes and wave functions of the pion and kaon with BHL prescription are proposed as the 
solution to the end-point singularities. The results show that the contributions from twist-3 parts 
are actually not power suppressed comparing with the leading-twist contribution. After including 
the effects from the transverse momentum of light meson valence-quark state and Sudakov factors, 
the decay widths of the Xcj into pions or kaons are comparable with the their experimental data. 
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I. INTRODUCTION 

The factorization form in the framework of the hard-scattering picture jl| is usually used 
in hadronic processes with large momentum transfer. In this picture, the full amplitude is 
factorized as a convolution of process independent distribution amplitudes of hadrons and 
process-dependent hard-scattering amplitude. However, the applicability of this approach at 
experimentally accessible momentum transfers, typically a few GeV, is questionable 2(]. One 
of the reasons is the large contributions from the soft end-point regions. A modified pertur- 
bative approach, so-call modified hard-scattering approach(mHSA), has been proposed by Li 
and Sterman where quark transverse momenta and Sudakov suppressions are taken into 
account. The advantage of this modified perturbative approach is the strong suppression of 
the soft end-point regions where the pQCD can not be applied. 

The exclusive charmonium decays have attracted interest for decades as they are an 
excellent laboratory for studying quark-gluon dynamics at relatively low energies. In the 
decay of P-wave charmonium Xco,2 to a pair of pseudoscalars, one finds that the lowest Fock 
state, the color-singlet contribution, alone is not sufficient to accommodate the experimental 
data. This discrepancy provides an important arena in which to test our understanding of 
the boundary domain between perturbative and non-perturbative QCD. From a naive point 
of view, QCD radiative correction is suppressed by the factor asiQ"^) and the contribution 
from higher Fock states, higher twist distribution amplitudes are suppressed by the factor 
l/Q"^- But this is not always the real case. In Ref. j^, the author showed that in the decay 
of Xco,2 PP (where P represents a light pseudoscalar meson.) the color-octet contribution 
from the higher Fock state contributes the same level as the color singlet state. 

To reveal the decay mechanism more clear, the systemic reanalysis on the contributions 
from higher twist distribution amplitudes will be significative and interested. When the 
amplitude of a physical process with large momentum transfer is only related to one 
hadron wave function, there is a suppression for the contribution from higher Fock states 
and higher twist distribution amplitudes. For example, the dominating contribution of 71 — 7 
transition form factor [sj] comes from the leading twist distribution amplitude of valence 
quark state. The QCD correction is only about 10% — 20% ^ and corrections from higher 
Fock state and higher twist distribution amplitudes Q] are suppressed by additional powers 
of l/Q^. However, for the exclusive processes with the overlap of the wave functions of 
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the initial hadron and final hadrons there are a lot of space to discuss the contributions 
from higher twist terms. For instance, the contribution of twist-3 distribution amplitudes 
to the pion electromagnetic form factor is comparable and even larger than the contribution 
from the leading-twist distribution amplitude of the pion at intermediate energy region of 
Q^, being 2 — 40 GeV^ in Ref. [8|]. The similar result has also been obtained for the 
kaon electromagnetic form factor in Ref. P]. More discussion about the contributions from 
the higher twist distribution amplitudes can be found in the Heavy-to-light transition form 
factors [l^, the nonleptonic two-body decays of the B meson 11 1, and so on. 

In this paper, we apply both the mHSA and the sHSA to reexamine the decay of 
Xcj n^TT^ , K^K^ (J = 0,2) including the contribution from the two-particle twist-3 
wave functions and distribution amplitudes of the light pseudoscalar meson. In the sHSA 
case, the end-point singularity can be avoided by using the BHL prescription [l^ in which 
the light meson distribution amplitudes are rewritten with exponential suppression factors. 
Our results show that the contributions from twist-3 wave function are comparable with or 
even larger than the leading twist contribution. Comparing with ^, the present analysis has 
an advantage that the theoretical uncertainty from light meson higher twist wave functions 
is less than that from the color octet wave function and the constituent gluon from the octet 
state of the charmonium. One can find the discussion on light meson twist-3 wave functions 

The paper is organized as following. In section II the main calculation of the hard- 
scattering amplitude in the modified perturbative QCD approach is presented. In section 
III, we present our model for the light meson two-particle twist-2 and twist-3 wave functions 
and the distribution amplitude within the BHL scheme. The section IV is the numerical 
analysis and section IV is our conclusion. The coefficients of hard-scattering amplitudes in 
the momentum space and in the b space are given in Appendix A and B. The Sudakov factor 
is presented in Appendix C. 



II. CALCULATION OF HARD-SCATTERING AMPLITUDES 

The two-meson decays of Xcj can be described as cc quarks annihilate into two gluons 
and then materializing into two final-state mesons as illustrated in Fig.l. We work in 
the rest frame of Xcj meson and take approximation M = 2mc where M is the mass of 
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(a) (b) 
FIG. 1: The hard scattering diagrams for the decay XcJ Tr'^Tr~ , K~ (3=0,2). 



the charmonium meson and rric is the mass of c-quark. The masses of pion and kaon are 
canceled in the Chiral hmit. Under these conventions, the momentums of the initial- and 
final-state mesons are described in terms of light-cone variables as 

, M M , M , , M , 

^^V^'V^'^^^' Pi = (^,0,Ot), P2 = (0,-^,0t), (1) 

where P is the momentum of the charmonium and pi and p2 are the momentums of 'k~^{K~^) 
and 7r~{K~), respectively. In the following, we only present the two-pion decay of XcJ for 
the two-kaon decay is similar. 

For the two-pion decay of XcJ, its decay width can be written as: 

Jz 

where A4 is the transitional matrix element. As usual, this matrix element within the 
mHSA can be factorized as the convolution with respect to the momentum fractions x, y 
and transverse separation scales bi, b2 of the two pions. 



MixcJ TT+TT ) = dxdy j 



A-n 

here the scripts = 7r,p, cr) mean the twist-2 and twist-3 wave functions for the pion 
meson and S is the Sudakov factor. T^j{x, y,hi,h2) is the Fourier transform of the hard- 
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scattering amplitudes T^j{x, y, ki, 

/^2l, ^2k. 

T^j^j{x, y, ki, k2) can be calculated from the graphs shown in Fig.l and the scripts of nonzero 
terms are ij = 7nT,pp,pa,(7p,cra, ki and k2 are the intrinsic transverse momentum of two 
final-state pion mesons, respectively. 

To calculate the transitional matrix element, the wave functions of the pion meson and 
the charmonium should be introduced. Similar to the definition of distribution amplitudes 
with leading and next-to-leading twist [IS], the light-cone wave functions of pion are defined 
in terms of bilocal operator matrix element 

< 71{P^)MZMZ2)\0 >= x/'^^'^/l^ g.(xp..., + (l-.K..,-k,.(z,-z,),) 

X {f^l5^n{x, kx) - /i,75(^p(x, kx) - CT^,p^;^{Zi - ^2) " ^"^g '^^^ ) ,(5) 

where /^r is the decay constant of pion and the parameter fi^, = ml / {rriu + mj) for charged 
pion. \l/7r, and are the twist-2 and twist-3 wave functions, respectively. The twist-3 
wave functions contributes power corrections. But at energy scale, the Chirally en- 
hanced parameter r^^ = f^n/fnc ~ 1 is large enough to consider the twist-3 contributions in 
charmonium decays. 

For simplif, the wave function of the charmonium is taken as [l^, 

X(P, q; J, J.) = J2 ^^'^^'^^ " < ^^'^'^'^^ > i^LM{q)VssAP, Q)^ (6) 

where q is the relative momentum between the quark and anti-quark, ipLMio) and VssXP^ l) 
are the spatial wave function and spin projection operators, respectively. The spin projection 
operators Vss^ {P, q) are 



'Pss.{P,<i) = y + (| + + ^)n55.K - (| - ^)] (7) 

with Uss^ = —75 for S* = and Uss^ = —^{Sz) for 5* = 1. Here s{Sz) refers to the spin 
part of the wave function, e{Sz) = (0,0,0,1) for 3^ = and e{Sz) = (0, ^1, -i, 0)772 for 
Sz = ±1. 

With the definition of the wave functions for the initial- and final-state hadrons, the 
transition matrix element of the two-pion decay for Xcj can be calculated in the coordinate 



space by standard method and then hard-scattering amphtudes T^j(x, ki, can be 
extracted as 

where li = xpi+yp2 + ^i—\<-2 and I2 = {l—x)pi + {l—y)p2 — ^i + ^2 are the momentum of two 
intermediate gluons, respectively. The operator 0^'^{q), which is related to the two-gluon 
annihilation of charmonium, reads 

where = {h — h)/'^ + q is the corresponding momentum of the c-quark propagator. The 
operators T^-J, are related to the two-meson materialization of two gluons with different twist 
wave functions and are expressed as 

= ^Tr [7^^^757.1*175], = :^^Tr [7,757.75], 



T---.f5^T.[7,757.a.,75].?.,^^^ 

Pa^ d d 



where the partial (^^ comes from the {zi — Z2Y term of Eq.(5). Here we don't take 

the momentum projection of Eq.(5), which can be obtained by transforming the parameters 
in terms of coordinate variable in Eq.(5) into the momentum space configuration, in Ref. 17 1 



or 



18|. 

For P-wave charmonium decay, the dominant contribution is given by q^ term of 



Tr[x(P, g; J, J z)0^'^ {q)] and the result after the integration of momentum q can be rewritten 

as 



'111 



E \/^^ {lM,lS-^JJ-^)e\M)Tr [of (0)^15. (P, 0) + a^'^(0)^i5.,A(P, 0) 
aT^Vvt 4 L 

with the definition of the partial of the P-wave function at the origin 

''^.V.m(.1 - ,/?^.\M), (12) 
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TABLE I: The values of coefficients C^, tj, dj, nj, ij in the eq.(14). 



I 



vrvr 



PP 



pa 



ap 



era 



ti 
di 
ni 
ii 



I 
4 

3 
1 
2 
1 



1 

3 

3 
2 
3 
3 



1 

3 

3 
2 
3 
3 



1 

9 

3 
3 
4 
6 



where e^{M) refers to the orbital part of the wave function and 



ar(o) 



(Zi-/2)V4-m2 ' 
rixY ih - k)xr[Ui - /2)/2 + m,]Y 



ih - /2) 74 - ml [ih - /2) 74 - mlf 



-P^sAm = W^K + |)(i + ^)ni5.K - f ) 



v.s.AP, 0) = J^[7a(i + ^)ni5.K - 1) + K + |)(i + ^)ni5.7A]. (is) 

The final hard-scattering amplitudes with transverse momentum can be expressed as 

TLMM) - ^,o-i^r-'''.''^"i0^±oHm'. (14) 

L- 12 i=Q 

where A^o = 167r'^/^/^mc^^-Rp(0), ao = 1 for Xco and a2 = for Xc2- The coefficients 

n/, i/ are listed in Table. I. The virtualities of the internal c-quark and the two 
intermediate gluons are 

Di = Amlxy — + ie, 

D2 = 4m^(l - x)(l - - + ie 

N = 2ml{x + y-2xy) + K^ -ie (15) 

with K = ki - ka. The coefficients CJ^(J) = 1 and Q^(J) = {-2y/^ml{x - yf are in 
agreement with the results from Ref. ^ and the others coefficients are showed in Appendix 
A. 
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Next, the fact that hard-scattering amphtudes T^j depend on ki and k2 only in the 
combination K imphes the following result for the Fourier transform of them 



dj— 1 dj—l nj—1 

1=0 1=0 i=0 

where jVo = l&T^'"^flB:^{f))/mT, = y^b, = ^y(l-x){l-y)b and rg = 
{x + y — 2xy)/2b with b = 2mcbi. H^p and Ki denote Hankel and modified Bessel func- 
tions, respectively. The (5-function, which simplifies the numerical work enormously, means 
that the two pions emerge from the decay with identical transverse separations. The coeffi- 
cients A{j(J), A2j(J) and B({J) are listed in Appendix B. 

The novel ingredient of the mHSA is the Sudakov factor , which takes into account 
those gluonic radiative corrections not accounted for in the QCD evolution of the wave 
function. In next-to-leading-log approximation, the Sudakov exponent reads 

S{x, y, bi, b2, h, ti) = s{x, hi, 2mc) + s{l - x,bi, 2mc) + s{y, 62, 2mc) + s(l - y, 62, SmJ 

1 , lii(fi/Ag,.^)hi(f,/Ag,.^) 

/3i ln(l/(6iAQCi.))ln(l/(62AQCD))' ^ ^ 

where the function s{x, bi, Q) with next-to-leading-log correction is given in Appendix C. 
The last term in Eq.(17) arises from a renormalization group transformation from the fac- 
torization scales fXpi to the renormalization scales tj at which the hard-scattering amplitudes 
T^j{x, y, bi, b2) are evaluated. 

The renormalization scales appearing in ctg and in the Sudakov exponent are chosen as 

ti = max{2mc^/xy, 1/bi, I/62}, 

t2 = max{2m^y/{l-x){l-y), 1/62} (18) 

by the virtualities of the intermediate gluons, which depend non-trivially on the integration 
variables. This choice of the renormalization scale avoids large logs from higher-order pQCD. 
The factorization scale is given by the quark-antiquark separation bj, fj,Fi = 1/bi. The ratio 
1/bi marks the interface between non-perturbatively soft momenta, which are implicitly 
accounted for in the pion wave functions, and the contributions from semi-hard gluons, 
incorporated in a perturbative way in the Sudakov factor. 



,2 
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Replacing e by 1 and ignoring the transverse momenta in T^j (x, kj, k2), one finds 
the decay amphtude within the sHSA as derived by Duncan and Mueller [l9|, 



A^(XcJ ^ TT+TT") = V / dx dy(j)i{x,fiF)T}jj{x,y,ml)(l)j{y,iXF), (19) 

Jo Jo 

where the renormalization scale tj is taken as the charm quark mass and customarily identi- 
fied with the factorization scale. The hard-scattering amplitudes T^j(x, y, ml) are expressed 
as follows 



T-r^ „ ^ ATpagK) x^ + {2- 6y)x + y{y + 2) 
iHo{^,y,mj ^^^^^ a;(l-a;Ml-y)(a; + y-2xy)2^ 

TMx,y,m^ 



-.pp t 2. _ A/'o"'("^c) f^l x2 + (3 - 8y)x + y{y + 3) 



72V6m2 ml x{l - x)y{l - y){x + y - 2xy) 

THo{x,y,m 



2, _ ^foal{ml) fil + x{-2y^ - 5y + 5)y + x\8y^ - 8y + 1) 



432\/6m2 x{l - x)y^{l - yy{x + y - 2xy) 

jfQ{x,y,m 



.ap , 2. _ ^oo^liK) A (1 - '^yW + y(% - 5)x2 + (5 - Sy)yx + 



JhoI^,Z/,^J - 2h92VQmlmlx\l-xYy\l-yy{x + y-2xyY^^'^y 

+1% - l)x^ + (48?/^ - 156?/3 + 1452/2 - 42?/ + 3)x^ + 2/(8 - llQy^ 

+218?/2 - 1292/ + 17)x3 + 2/(-122/^ + 892/^ - 972/^ + 242/ + 4) x^ 

V(62/' - 262/2 + 92/ + 4)x - (2/ - 3)2/') (20) 



for J = and 



T-rr 7, 7772^ - ^0"'^"^') X + 2/ - x2 - 2/2 



36y3m2 x(l - x)2/(l - y){x + y - 2xyf 

T^P^i^^y^ml) 



2 A/'oa2(m2) ^2 _ 



72V3m2 ml x(l - x)2/(l - 2/) (a; + '2/ - 2xy) 

T^2(a;,l/,"^c) 



2, _ Aroa2(rri2) ^2 ^§^2 _ ^ ^^^2 + 22/(-42/2 + 22/ + l)a; + 2/2(42/ - 3) 



432\/3m2 ml x(l - x)2/2(l - yf{x + 2/ - 2x2/)^ 
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Moal{ml) (4 - 8y)x^ + {Sy"^ + ^y - S)^^ + 2(1 - Ay)yx + y'^) 



432y3m2 ml 



x'^il — xYyil — y){x + y — 2xyY 



TH2{x,y,ml 



for J = 2. 



-((24|/3-36y2 



2592y3m2 x^{l - x^y^il - yy{x + y- 2xyY 
+lAy - l)x^ + (-144y^ + 228y^ - llV + 30?/ - 3)a;^ + (8 + 268?/^ 

-424y^ + 210?/^ - 46y + 3)x^ + y{-l2y^ - 175y^ + 2A2y^ - 84 y + 10)a;^ 



+2|/2(3|/3 + 23y^ - 27y + 5)x - /^{y^ + 3y - 3)) 



(21) 



III. THE WAVE FUNCTIONS OF LIGHT MESONS 



In the above calculation, the twist-2,3 wave functions and distribution amplitudes of pion 
and kaon are the main non-perturbative input parameters for mHSA and sHSA, respectively. 
In this section, we will discuss them in detail. According to BHL prescription 12j, one can 



connect the equal-time wave function in the rest frame and the light-cone wave function by 
equating the off-shell propagator in the two frames. The wave function for quark- ant iquark 
systems at the infinite momentum frame can be got from the harmonic oscillator model at 
the rest frame 



\E'(x, k_|_) oc exp 



k^ + mj k^ + ml 



X 



+ 



1 — X 



(22) 



where mj and (3 are the constitute quark mass and the harmonic parameter, respectively. 
The distribution amplitude can be obtained from the integration of wave function over the 
transverse momentum 

rf^k. 



(X^flp) 



IGtt' 



-^(a;,kx) 



(23) 



\kj_\<fj.p 

where /Xi? is the ultraviolet cutoff. 

Upon expansion over Gegenbauer polynomials, twist-2 wave functions of pion and kaon 
with the transverse momentum dependence can be characterized as 



^^(x, kx) = Al 1 + BlCl'\2x - 1) + ClCl''\2x - 1 



-y3/2. 



exp 



\^\ + ml 



W^xil 



X] 



(24) 
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and 



xexp 



+ 



X 



X 



(25) 



where Cn^^ arc Gcgcnbaucr polynomials and q means light quark u or d. For the SU(2) iso- 
topic symmetry, the odd expansion terms do not appear in the pion wave functions. On the 
contrary, the odd expansion terms are not zero in the Kaon wave functions for SU(3)-fiavor 
symmetry breaking. Estimates of first two Gegenbauer moments for twist-3 distribution am- 
plitudes are more uncertain than that of leading twist distribution amplitudes. To simplify 
the following numerical analysis, we take twist-3 wave functions as 



x(l — X 



-exp 



Plx{l-x) 



, ^l{x,k^) = Alexp 



- x) 



(26) 



and 



x{l — X 
^^{x,k^)^A^exp 



rrig kl + mj 
1-x 




(27) 



for pion and kaon, respectively. 

Substituting Eq.(24)-Eq.(27) into Eq.(23), the distribution amplitudes of pion and kaon 
are written as 



- x) [l + BlcT{2x - 1) + C:cT{2x - 1)] exp 



SI3lx{l-x) 



,(28) 



^k{x) = 4f^^(l -x)\l + Bfcl'\2x - 1) + C^cT{2x - 1) 



27r2 

xexp 



[1 — x)m'i -\- xm 



Wkx{1 - x) 



for twist- 2 distribution amplitudes and 



(t>l{x) = -^exp 



3^x{l — x) 



Kg 

27r2 



■exp 



3^x{l — x) 



(29) 



(30) 



K (^) = ^^^^^ 



1^1 — x)mg + xml 



3lx{l 



x) 



27r2) 



-exp 



[1 — x)rn?q + xml 
^I3lx{l - x) 



(31) 
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for twist-3 distribution amplitudes. With the help of the above distribution amplitudes 
from BHL prescription, the endpoint problem can be cured in the standard HSA since the 
exponential suppression appears in x = and x = 1 point. 

For definiteness, we take the conventional values for the constitute quark masses: rrig = 
O.SOGeV and rrig = OAdGeV. The parameters, A^, Bj, Cj and Pj {i = vr, K,p, a; j = vr, K) 
can be determined by some constraints on the general properties of the light mesons wave 
functions. In the pion and kaon case, the harmonic parameters /J^r and Pk are obtained 
by the constraints < >j^^< >7r~ (0.356GeV^)^, which are the average values of the 
transverse momentum square defined as 

< kl >,= I J dx^|kl||vl/,(x,k^)|VP;- (32) 

with i = TT,K and stand for twist-2 wave functions \E'^ and The decay constants are 
taken as /tt = 0.132 GeV for the pion and fx = 0.160 GeV for the kaon. The probability of 
finding the qq leading-twist Fock state in a pion or kaon should be not larger than unity, 

= ^ / dx^\^,ix,k^)\' < 1. (33) 



24 J 167r3 

The others coefficients are extracted by the normalization condition 

dx^^l{x,k^) = 1 (34) 

with [i = TT,K,p,a] j = tc,K) and first two Gegenbauer moments of twist-2 distribution 
amplitudes aj (z = 2, 4 for j = tt; z = 1, 2 for j = K.). The coefficients at some reference 
scale lip are nonperturbative quantities and have to be evaluated using a nonperturbative 
technique or must be extracted from experiment. It turns out that the determination of a 2 
receive large errors, whether by direct calculations using QCD sum rulesj20] or by analysis 
of experimental data on the pion electromagnetic and transition form factors 2l|| . Totally, 
the averages of the second moment are probably 

4(1 Gel^) = 0.25 ±0.15, af(l GeV) = 0.25 ± 0.15, (35) 

in Ref . |l^ , including radiative corrections to the sum rules. 

The numerical value of the first moment af was the subject of significant controversy 



until rece ntly . 



rules 



22 



23 



24 



le existing estimates are all obtained using different versions of QCD sum 



251 ] and yield an average [li 

af (1 GeV) = 0.06 ± 0.03. (36) 
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TABLE II: The parameters of twist-2,3 wave functions for the pion and kaon mesons in Eq.(24)- 
Eq.(31). The dimensions of harmonic parameters (i-w^K and normahzation coefficients (i = 
Tr,K,p,a;j = tt^K.) are GeV and {GeV)~'^, respectively. The others parameters is dimensionless. 









TT 


TT 


TT 


p 


a 


upper hmit 


0.512 


0.270 


672.28 


0.628 


0.354 


106.92 


574.47 


central value 


0.461 


0.249 


849.18 


0.469 


0.213 


140.88 


747.55 


lower limit 


0.418 


0.259 


1034.30 


0.259 


0.023 


184.29 


965.85 


K 


Pk 


pK 

m 








Zip 




upper limit 


0.461 


0.492 


1108.83 


0.253 


0.618 


170.59 


883.49 


central value 


0.442 


0.422 


1196.66 


0.215 


0.477 


193.81 


998.32 


lower limit 


0.417 


0.398 


1353.56 


0.175 


0.323 


232.44 


1188.23 



Estimates of yet higher-order Gegenbauer moments are rather uncertain. The fourth Gegen- 
bauer moment of the pion twist-2 distribution amphtude [261] was constrained 

<(1 Gel/) = 0.04 ±0.11, (37) 

which is consistent with the results from the hght-cone sum rule calculations of the transition 
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28 



29| 



form factor F^^^^* in Refs. 

According to QCD evolution of the wave function, the coefficients a] at a factoriza- 
tion scale fip can be expressed as a^fip) = al{fil){as{fJ''j?) /as{f^o)V' ■ cti(/^o) ^ non- 
perturbative coefficients, /xq is a typical hadronic scale, 0.5 < fio < I GeV, and 7, are the 
anomalous dimensions. In this work, a reasonable factorization scale should be chosen as 
fip = rric, 1.35 < rric < 1.8 GeV. In Ref . |l^ , the values of the coefficients a{ at /i^ = 1 
GeV and /ii? = 2 GeV are listed in Table. 3. For example, ai{lGeV) = 0.06 ± 0.03 and 
ai{2GeV) = 0.05 ± 0.02. By analyzing these data, we find that it is feasible to choose 
Gegenbauer moments at fip = 1 GeV in our calculation. 

Taking account of the above Gegenbauer moments for the pion and kaon twist-2 distri- 
bution amplitudes, we figure out harmonic parameters Pn,K, probabilities of finding the qq 
leading- twist Fock state P^^^ and Gegenbauer coefficients Ai, Bj, Cj {i = n,K,p,a; j = 
7r,i^). Those values are list in Table. II. According to uncertainties of twist-2 Gegenbauer 
moments, the parameters of the pion and kaon are given in three parts: upper limit, central 
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value and lower limit. Since P^^ < 0.270 and P^^ < 0.492 are much smaller than unity, 
higher twist and higher Fock states are important components of the pion and kaon. 

In the mHSA, the convolutions of wave functions and hard- scattering amplitudes are 
presented in transverse configuration b-space. We need to define wave function in b-space 
by Fourier transformation 

m\x, h)= [ k^)e-^'^^-^ = 47r0*(x)exp \-2Pfx{l - x)b'^] , (38) 

J (27r)^ 

where and 0* stand for wave function and distribution amplitude 0* {i = vr, K, p, a; j = 
IT, K), respectively. One may observe that wave functions in the b-space are also highly sup- 
pressed in the endpoint region. Such feature is necessary to suppress the endpoint singularity 
coming from the hard-scattering amplitudes and then to derive a more reasonable results. 



IV. NUMERICAL ANALYSIS 



In our calculations for the decay ratio of Xcj to light pseudoscalar pairs, the partial of 
P-wave function at the origin -Rp(O) is also an important nonperturbative input parameter. 
It is shown that this parameter is a function of the charm-quark mass both in the well- 



cnow quarkonium potential models 



30 



3l| and in the global fit of charmonium parameters 



32l | . To obtain its expression relative to the charm-quark mass, which is consisted with our 
approach, the decay width of the Xco annihilating into two photon need to be calculated by 



the same approach. With the help of Refs. 
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33 



34 



4„ 2 

em 



\r: 



mi 



35| . we obtain 

2 28.aJm2l 



TT 



(39) 



I is the charge of the c-quark 



where the one-loop QCD radiative correction is included, Qc 
and = 1I7 is the electromagnetic coupling constant. 

The running coupling constant asiQ"^) up to next-to-leading-log is written as 

asm 1 P,\tiHQ^/KIcd) 



TT 



(40) 



- PMQyAl^^) Pf ln^(QVA|cD) 
with f3i = (33 — 2n/)/12 and /?2 = (153 — 19?2/)/24. Here we take quark-flavor number 
Uf = 4 and the QCD scale Aqcd = 0.25 GeV. 



The decay width of the Xco annihilating into two photon can be obtained from Refs. 36|, 



371]. Using the above data and formulas, the relation between -R' (0) and rric is shown in 
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Fig. 2 witli uncertainty of tlie Xco two-plioton decay widtli. Tlie solid curve comes from 
taking central value T[xco 77] = 2.37 keV. The dot-dashed and dashed curves are given 
by taking upper limit T[xco — ^ 77] = 2.71 keV and lower limit T[xco — ^ 77] = 2.03 keV, 
respectively. We find that uncertainty with different Xco two-photon decay width is less 
than 10%. So we will take the result of the central value in the following. The region of the 
charm-quark mass is rric = 1.35 — 1.8 GeV from Ref.[J]. Comparing our -R1(0) with values 
of Ref. ^, there are some differences that our value is less than one of Ref. ^ as rric = 1.35 
GeV and vice versa as rrir =1.8 GeV. 



> 




nic GeV 

FIG. 2: Dependence of the partial of P-wave function at the origin i?p(0) for P-wave charmonium 
on the c-quark mass rric with uncertainty of the XcO two-photon decay width. 

On the other hand, the chiral enhancing scales /^tt and which are scales characterized 
by the chiral perturbation theory, are important parameters which can affect contributions 
from twist-3 parts sensitively. However, they are difficult to give precise numbers as long as 
the current quark masses are not more accurately known. To obtain reasonable numerical 
analysis with acceptable estimates for the chiral enhancing scales, we take /i7r(l GeV) = 
1.5 GeV andu^l GeV) = 1.7 GeV, which are consistent with the results from pQCD 



application 
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391 ] and chiral perturbation theory [40]. 
In the above discussion, nonperturbative input parameters appearing in our calculation 
are obtained by a model wave function or distribution amplitude, a fit experimental data 
and a reasonable evaluation. Next, we will do numerical analysis for the charmonium Xco,2 
decay rates into two pions or two kaons with the charge case in the sHSA and in the mHSA, 
respectively. 

The decay widths of XcJ tt+tt^ and K^K^ (J = 0, 2) by the sHSA are shown in Fig. 3 
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1.4 1.5 1.6 1.7 1.8 1.4 1.5 1.6 1.7 1.8 



iMc GeV trie GeV 




1.4 1.5 1.6 1.7 1.8 1.4 1.5 1.6 1.7 1.8 



Me GeV trie GeV 

FIG. 3: Dependence of the prediction for the XcO ~^ vr+vr^ and K^K~ decay widths on the c- 
quark mass rric with contribution from twist-3 distribution amphtudes of the pion and kaon meson 
in sHSA, respectively. 



and Fig. 4 with the c-quark mass rric as a variable parameter. We take the central-value of 
Table. II as the input parameters for the top two and middle two figures. The top two figures 
are the decay amplitudes for two-pion and two-kaon decays, where the dashed curve is the 
contribution from twist-2 part, the dotted and dot-dashed are contributions from twist-3 
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FIG. 4: Dependence of the prediction for the Xc2 vr+vr^ and K^K" decay widths on the c- 
quark mass rric with contribution from twist-3 distribution amphtudes of the pion and kaon meson 
in sHSA, respectively. 



parts named a and b and the solid curve is the sum. The twist-3a is the positive and the 
twistSb is the negative which are corresponding to the positive and negative terms in Eq.(20) 
and Eq.(21). 

The large enhancement of the total amplitudes in the X2 channel indicates that the twist- 
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TABLE III: The decay widths for XcJ ~^ vr+vr^ and K^K^ (J = 0,2) from experimental data. 
The BES results are evaluated with the BES result for the total width. In the other cases the PDG 
average values for the total widths are used. 









PDGf33 


BES[4]J 


Belle [42] 


r[xco - 




TT ] [kev] 


50 ±8 


67 ±36 


60± 21 


r[xco - 


^K^ 


'K-][kev] 


60± 10 


81 ±45 


57± 19 


nxc2 - 


-^7^^ 


"vr^] [kev] 


2.8 ±0.5 


3.0 ± 1.0 


3.4 ± 1.3 


nxc2 - 




^K-][kev] 


1.5 ±0.4 


1.6 ±0.7 


2 ±0.9 



3 distribution amplitudes play important pole. The results for the decay width are shown in 
the middle two figures of Fig. 4 where the solid curve is the total decay width and the dashed 
curve is the for twist-2 part. In the Xo channel, the corrections are not so large since the two 
contributions from twist-3a and b parts have opposite sign and have large cancelation. The 
results with the uncertainty of twist-2 distribution amplitude are shown in the bottom two 
figures, where the solid curve is for the central value, the dashed and dot-dashed curve are 
for lower limit and upper limit, respectively. We can see the uncertain is very large which 
shows the sensitivity on the distribution amplitude. 



Comparing the results with experimental data 37|, |4l|, |42l] , which are list in Table. Ill, we 



see that the decay widths of Xo channel are smaller than experimental data for all variable 
rric and the decay widths of the X2 channel are in agreement with the experimental data in 
the region rric < 1.5 GeV. 

In Fig. 5, we show the curves of the widths of the Xco,2 into two pions or two kaons on 
the c-quark mass mc by the mHSA method. The solid curve is the decay width where only 
the twist-2 contribution is considered and the parameters are taken as the central values. 
The shadow is the total decay width where the uncertainty of twist-2 Gegenbauer moments 
are considered. Here we see that the shade region is very narrow which means the un- 
sensitivity on the twist-2 distribution amplitude of the light mesons. The decay widthes 
including twist-3 corrections are improved remarkably. Detailedly, the predictions for the 
decay widths of X2 tt+tt^ and xo K^K~ are comparable with experimental data in 
the region ttlc G (1.4, 1.6) GeV and rric G (1.35, 1.6) GeV, respectively. This is very different 
with the sHSA method which suggests the necessary of the mHSA method. 
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V. CONCLUSION 



In this paper, we presented a detailed analysis of Xco,2 decays into two pions and two 
kaons including the twist-3 contribution within the framework of the sHSA and mHSA 
methods. In the sHSA, the end-point problem is overcame by using BHL prescription where 
a exponential suppression is introduced in the expression of hadronic wave functions or 
distribution amplitudes. The uncertainty of the results on the twist-2 Gegenbauer moments 
for the pion and kaon is analyzed and is rather small in the mHSA method. The results 
indicate the larger contributions from twist-3 distribution amplitude which have not been 
analyzed before. And both the decay widthes of Xco,2 to tt+tt^ and K^K~ are found to be 
comparable with the experimental data in the region nic € (1.35, 1.8) GeV when including 
twist-3 correction in the mHSA. 





1.4 1.5 1.6 1.7 1.8 1.4 1.5 1.6 1.7 1.8 



trie GeV iHc GeV 

FIG. 5: Dependence of the prediction for the XcJ vr^vr^ and K^K^ (J=0,2) decay widths on 
the c-quark mass rric with contribution from twist-3 distribution amplitudes of the pion and kaon 
meson in mHSA, respectively. 
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APPENDIX A: THE COEFFICIENTS OF HARD-SCATTERING AMPLITUDES 



In this Appendix we present the exphcit expression of the coefficients C/(J)(/ = 
pp,pcr, ap, aa) appearing in the Eq.(14) for hard-scattering amphtudes. 

C7f(0) = 5, C^^{0)^2ml{x' + 2{y-l)x + {y-2)y), (Al) 
Cf (2) = 2, C^(2) = -Amlix^ -Ayx + x + y^ + y) (A2) 

for hard-scattering amphtude T^j{x, y, ki, with J = 0, 2. 

Cf (0) = bml Cr(0) = 4m^(4x^ + {Ay - 6)x - 2y - 3), 
Cf''(0) = -Amlillx^ + 2(7?/ - 9)x + ?/(3?/ - 10)), 

CP'^(O) = -16m^(2a;^ + {<6y - 7)x^ + 3{2y^ - 5y + 2)x^ + y(2y2 - 9y + 8)x 

- (y - 2)y^) (A3) 

and 

Cf (2) = 8m^, Cr(2) = -8m^(4x' - 8yx + 4^/ + 3), 



Cr(2) = 16m^(4x' + (3 - 14y)x + 7y), 

Co^"(2) = 64m^(x^ - 2x' - ^y - l)yx' - 2y{y'' -3y + l)x + {y - 2)y^) (A4) 

for hard-scattering amphtude T^j{x, y, ki, k2) with J = 0, 2. 

Cl^Qi) = hml Cr(0) = 4m^(4y2 - 6y + x(4y - 2) - 3), 

C7(0) = -4m^(3a;2 + 2(7?/ - 5)x + |/(11|/ - 18)), 

C^P(O) = -16m^((2y - l)^^ + (6y2 - 9?/ + 2)x^ + y(6y2 - 15y + 8)x 

+ yW-7y + 6)) (A5) 

and 

Cr{2) = 8ml, Cr{2) = -8ml{Ay' + x{A - Sy) + 3), 
Cr(2) = 16m^(x(7 - 14y) + y(4y + 3)), 

C^P{2) = -64ml{{2y - l)x' + {3y^ - 6y + 2)x^ + (2 - 3y)yx - {y - 2)y^) (A6) 
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for hard-scattering amplitude T^j{x, y, ki, with J = 0, 2. 

Q<^(0) = 7ml, Q"(0) = -2m^(7a;2 -lx+ lly^ - lly + 39), 

CT(0) = 2m^((26 - 52y)x^ + (-lOV + 182y - 101)x^ - 2(26y^ - 91y^ + 173 y - 85)x 
+ 26y^-6V + 138y + 40), 

C^'^(O) = -16ml{x^ + (8y - 6)x^ + 2{5y^ - lly - 3)x'^ - 2{5y^ + ISy - 18)x - 3 

+ 67/3-6y2 + 28t/-6), 
C^^(O) = 32mJ°((8y - 4)x^ + 4(8y^ - 13y + 5)x^ + 2(24y^ - 68y^ + 62y - 19) 

+ {32y^ - 136y^ + 200y^ - ISOy + 35)x^ + 2{Ay^ - 26y^ + 5Ay^ - 53 y^ + 31y 

- 9)x + yi-Ay"" + 12y^ - Uy^ + 19y - 18)), 
Cf^(O) = 32mJ^(x^ + (1 - 8y)x^ + (-45^^ + 65y - 9)x^ - 2{A0y^ - lOSy^ + 50 y + 6)x^ 

+ {-65y^ + 250y^ - 2Uy^ - 36y + 36)x^ + y{-2Ay^ + 125y^ - 164 - 20y 

+ 72)a; + y^(-3y^ + 21^^ - 41y2 + 4y + 36)), 
Co"^(0) = -96mJ"(x2 + 2(y - l)x + (y - 2)yf{-2y + x{Ay - 2) - 1) (A7) 



and 



C^"(2) = Ami, ^""(2) = 4m^(7,r2 - (12y + l)x + lly=^ - 5y - 6), 

CI" {2) = Aml{26{2y - l)x^ + {-2Q8y^ + 130y - 61)x'^ + {52y^ + 130y^ - 32 y + 25)x 

-26y^-93y2 + 57y-4), 
CT(2) = 32m^(a;^ + (21 - 46y)a;3 + (118y2 - 49y + 15)x^ - y{bAy'^ + 37y + 6) a; - 3y^ 

+ 33y^ + 15y2-8y + 6), 
C^<^(2) = -32mJ°(8(2y - l)x^ + (-32y2 - 8y + 7)x^ + (-96y3 + 208y2 - 166 y + h3)x^ 

+ (-32y^ + 208y^ - 26y^ - 29y + 10)x^ + (16y^ - 8y^ - 198 y^ + 19y^ 

- 56y + 30)x + y(-8y^ - 9y^ + lOly^ - 22y + 30)), 

Cf'^(2) = -Umf{x^ + (13 - 32y)x^ + (Sly^ + 29y - 12)x^ + 2(80y=^ - 171 y^ + 64y 

- 12)a;^ + (31y^ - 302y=^ + 248y^ + 24y - 24)x^ + y(-48y^ + 89 y^ + 64y^ 
+ 40y - 48)a; - y^(3y^ - 33y^ + 44y^ + 8y + 24)), 

C^''{2) = 192ml\x^ + 2(y - l)x + (y - 2)y)2((4y - 2)x^ + (-16y2 + lOy - 3) + {Ay^ 
+ lOy^ - l)x - y(2y2 + 3y + 1)) (A8) 
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for hard-scattering amplitude T^j{x, y, ki, with J = 0, 2. 



APPENDIX B: THE COEFFICIENTS OF THE FOURIER TRANSFORM 

In this appendix we present the exphcit expression of the coefficients A{j(J), A{^{J) 
and Bj (J) appearing in the Eq.(16) for the Fourier transform of hard-scattering amphtude. 
First, we show some useful formulas for Fourier transform, 

d'\^ e~^^-^ _ (-r \ i£^[<\V~sb)] for.>0 ^^^^ 



(27r)2 (, - + ^e)n (n - 1)! | ^£^[MV^sb)] for 5 < 
with n = 1, 2, • • • . 

The coefficients A{-{J), A^.(J) and S/(J), coming from the Fourier transform of hard- 
scattering amplitudes, are given as, 

.^^.j._ 2{x + y) + {-2y/\x-yr 2(2 - x - y) + i-2y/\x - yf 

{x + y){2 - X - y) + {-2y/^{x - yf 



■k{x + yY{2 - X - yy 



4v27r(x + y){2 - X - y)\/x + y - 2xy 
with b — 2mcbi for T^j{x, y, bi, b2). 



A{x + y - l){x + yY ^ A{x + y - l)(2 - x - y)' 



B^(0) = - -, Sf(0) = , (B3) 

Ti{x + y){2 - X - yY 2^ix^x + y - 2xy ^ ' 

for r^^(x,y,bi,b2), 

.PP.ox _ 4a;y - - _ 4(1 - - y) - [l - xY - (1 - I/)' 

^^°^^^-2(x + y-l)(x + y)2' ^2°^^^- 2(x + y-l)(2-x-y)2 

7r(x + y)2(2-x-y)2' V27r(x + y)(2 - x - y) + y - 2xy ^ ' 
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for T^^{x,y,hiM), 

_ + {2y - 3)x -2y' + y ^ + {2y -7)x-2y' + y + 2 

^ 8ix + y-mx + yr ' ^ 8(x + y - 1)2(2 - x - y)3 ' 



16(x + y-l)(x + y) I6(x + y - 1)(2 - x - y) ^(1 - x)(l - y) 



^0^^°) - 32.(. + y)3(2 - . - y)3(. + y - 2.y) ^^^^^^ + ^^^^"^ - ^^^^^ + ^^^^^^^ " 
+4)x^ + 2((10y^ - 35 y^ + 32y - 6)6^ - 160y + 80)x=^ + 2(56^y^ - 25b'^y^ 

+4(96^ - 16) y^ + (304 - 146^)^ - 96)x^ + y(262y4 _ 1552^3 ^ 32(52 ^ g)^2 
-20 (6^ + 8)y - 128)a; - y^{b'^y^ - Ah'^y'^ + 4(6^ + 24)y - 64)) 

_ ^(2a^^ + (2% - 17)x3 + (14y2 - 53y + is)^^ + y(-10y2 + y + I6)x + y2(5 y - 2)) 
8v 27r(2 — a; — y)2(a; + yY\/x + y — 2xy{2xy — x — y) 

B^'^iO) - 6^(2x^ + (2y-3)x-y) 

^ ^ ^ 327r(a; + y)(2-a;-y)(a; + y-2a;y) ^ ^ 

for r|o(a;,y, bi,b2), 

,pa.2^ = 43:- + i2fj 3).r -2tj^ + y ^ Ax^ + {2ij 7).r -2tj^ + tj + 2 

^ 4(x + y-l)2(x + y)3 ' ; 4(^ + ^_ 1)2(2 -x-y)3 ' 

,^pa(2) - ^^^^ ~ ^2^(2)- 6(x-l)(x-2y + l) 



8(x + y - l)(x + y)2^' 8(x + y - 1)(2 - X - y) V(l - - y) 



^0^^^) - 8.(. + y)3(2 - . - y)3(. + y - 2.y) ' '"'^^ ' + ^^^^^^ ' 

+4(26^^3 _ 4^2 y2 ^ (^2 ^ 4Q)^ _ 20)2;3 ^ (752^4 _ 26^,2^3 ^ 3(3^2 ^ g)^2 

-4 (6^ + 76)y + 96)x^ + 2y(62y^ - 662y^ + 2(56^ - 24)y2 - 4(6^ - 10) y + 32)x 
V(&V-46V + 4(6^-12)y + 32)) 

^p.(2) = -Ka^' + 2(8y - 5)x3 + (-11^^ - 13y + 6)x-^ - 2y(y2 - 7y - l)x + {y - A)y^) 

2\/27r(2 — x — yY{x + yY^/x + y — 2xy{2xy — x — y) 

B^''i2) = -&'(x^-2xy + y) ^ 

^ ^ ^ 87r(a; + y)(2-a;-y)(a; + y-2a;y) ^ ' 

23 



for Tg'^(x,y,hiM), 

^ 2a:^ - {2y + l)x + (3 - 4y)y ^ -2x^ + 2yx + x + Ay^ - 7y + 2 

^ 8{x + y-mx + yr ' ^ 8(x + y - 1)2(2 - x - y)3 ' 



16(x + y-l)(x + y)^' I6(x + y - 1)(2 - x - y) ^(1 - x)(l - y) 



^0^^°) - 32.(. + y)3(2 - . - yfix + y - 2.y) ^^^^^^ " ^^^^ + ^^^^"^^ " + 

+2((10y^ - 25y^ + 16y - 2) 6^ + 96y - 48)x^ + 2{m^y^ - 35b'^y^ + 4(96^ 

-16)y^ - 10(6^ + 8) y + 32)x^ + y(106V - 456^?/^ + 64(6^ - 5)y2 + (608 

-2862)y - 128) x + y^{2b'^y^ - 116V + 206V - 4(36^ - 40)y - 192)) 

^ap^Q^ ^ 6((5 - Wy)x^ + (14y2 + y - 2)x^ + y(26y2 - 53y + lQ)x + y2(2y2 - 17 y + 18)) 

8\/27r(2 — a; — y)2(a; + yY\Jx + y — 2xy{2xy — x — y) 

^.WQ^ ^ 6^(y(2y-3)+a;(2y-l)) 
^ ^ 327r(a; + y)(2-a;-y)(a; + y-2a;y) ^ ' 

for 7^g(x,y,bi,b2), 

p ^ -2.7-^ + 2</.r + .r + </(!</ - 3) .ap/oN ^ -2.7-- + 2tj.v + .)■ + -ly + 2 
^ 4(x + y-l)2(x + y)3 ' 2oU + ^ _ 1)2(3 _ a; _ y)3 ' 

A^f (2) = ^(^ - 2a:) - - l)(y - 2x + 1) 



8(x + y-l)(x + y)2^' 8(x + y - 1)(2 - X - y) V(l - - v) 



= ^ , ^ ,3^9 --{h\2y-l)x' + h\ly'-l2y + A)x' 

87r(a; + y)'^(2 — x — yY[x + y — 2a:y) 

+((8y3 - 26y2 + 20y - 4)6^ - 96 y + 48)x^ + 2(62y^ - Sh'^y^ + 4(36^ + 8)y2 
-4(6^ - 10)y - 16)x^ + y (-262y^ + 36^y3 + 4(6^ + 40)y^ - 4(6^ + 76)y + 64)x 
-y^(62 y^ - 462y^ + Aly^y"^ + 80y - 96)) 

^ap(2) = K(2|/ - 1)3^'^ + (1 V - 14y + 4)a;=^ + y(-16y^ + 13y - 2)x - y\y^ - 10 y + 6)) 

2\/27r(2 — X — y)^{x + y)^\/x + y — 2xy{2xy — x — y) 

B''P(2) = -b\y^-2xy + x) 

^ ^ ' 8n{x + y){2 - X - y){x + y - 2xy) ^ ' 
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for 7;jf(x,y,bi,b2), 

^10 (0) = i28(x + -l)3(x + y)^ ^^'^' + ^^^^ " " + + - % - 1) 

-8(48?/^ - 36 y + + (b^(20y^ - lOy^ - % + 1) - 48(12y=^ - lOy^ + Ay 

+l))x^ + (3 (56^ - 128)i/^ + (352 - m'^)y^ - 2{3b'^ + 88)y^ + 3(6^ + A8)y 

+24) + y{6{b'^ - 16)t/^ + (96 - 5b'^)y^ - A{b'^ - A8)y^ + 3(6^ - I6)y - 80) x 

+y\bY - bV - {h' - + {b' - 240)y + 88)) 

^20^(0) = i28(x + y - 7)3(2 - x - y)^ ^^^^^^^ " + ^^^^^^ " + ^^^^^ + 

-1747/2 + 240y - 67) + (48^/^ - 3647/3 + 9347/2 _ 9267/ + 24Q)x^ + 2(67/^ 

-72y^ + 280 y^ - 499^^ + 397y - 100)x - I2y^ + 73y^ - 206^^ + 305^^ - 216y 

+52) -b'^ {x + y -3){x + y - 2f{x + y- if) 



(0) = ^T^— ^ .3 ^ ((1 - 87/)a;^ + (-247/^ _ 28y + 6)a;=' - {24y' + 587/^ 

64(x + 7/ - iy{x + yr^/xy 



-34y + 3)a;2 + 2y (-4y=^ - Uy^ + 17y + l)x + y^{y^ + 6y - 3)) 

^21^(0) = -— -— ^ .3... ... =((7 -87/)x^- 2(127/2 -547/ + 37)^3 

64(x + y — 1) (2 — a; — T/j^y/^l — a;)(l — y) 

+ (-247/3 + 2027/2 - 4307/ + 231)^2 - 2 (47/^ - 547/^ + 2 1 57/2 - 3077/ + 134)x 
+7y^ - 74y3 + 231y2 - 268y + 100) 



128(x + 7/)(x + 7/- 1)' ' 128(2 - X - y){x + y - I) 

^"'C" = 256.(. + .)^(2 - . - yn. + , - 2.yy- ""^^^^ ' + ^"''^'^'^^ " ^^^^ 
+121) - 128(1 - 2yf)x^ + 8 ((1987/^ - 4657/2 + 3097/ - 59)^2 + 16(-167/3 

+407/2 - 267/ + 5))a;^ + 4 (62(8927/^ - 3170^^ + 35717/2 - 14827/ + 186) 

-32(12y^ - 80y3 + 72 y2 _ _ + 4((1110y^ - 5451y^ + 90067/^ 

-60067/2 ^ ^r^^^y _ gg^ ^2 ^ 32(32^5 _ 44^^ _ ^-^^y^ ^ 45^2 ^ 34^ _ ^2))x^ 

+2((15607/^ - 10230 7/^ + 232357/^ - 222607/^ + 86207/2 - 96O7/ - 24)^2 
+64(68y^ - 284y5 + 246 y^ + 203^^ - Il5y2 _ 28y + 12))x^ + 8y((142y^ 
-1277y^ + 3947 - 5225^^ + 2868y2 - 464y - 24)^2 + 16(48y^ - 304^^ 
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+730y^ - 537 - 196y^ + 168y - 16))x^ + V((36/ - 570y^ + 2603y^ 
-4878 + 3782?/^ - 86% - 72)6^ + 32(12?/'' - 120?/^ + 496?/^ - 967?/^ 
+765 y2 - 88y - 24))a;2 - 2y3((6y6 + 45^5 _ Q^2y'' + 1924y3 - 2264y2 + 768 y 
+96)6^ + 64(12y^ - 78y^ + 281y^ - 482^^ + 348y - 48))x + y^(6^(6 y^ + 33?/2 
-76|/ - 12) (1/ - 2)2 + 128(3|/^ - 15?/^ + 5^ - 76|/ + 44))) 

5r(0) = ^ (-66^(1 -2y)V-(2|/-l) 

512V27r(a; + y)3(2-,T-?/)3(x + |/-2x|/)5/2^ 

(36^(24^2 - 40y + 11) - 16(32y2 - 32 y + 13))x^ + (8(512/ - m2y^ + 1224^^ 
-446y + 61) - 36^(120?/^ - 408 y'^ + 428?/^ - 161?/ + 18))x^ + (16(384y^ 
-1728|/^ + 33181/^ - 2541 y'^ + 863y - 103) - 36^(160?/^ - 760y^ + 1186i/^ 
-715t/2 + 148y - 4)) + (8(5127/6 - 3456?/^ + 10480?/^ - 165907/^ + llOlly^ 
-2982y + 236) - 3 h^{l2Qy^ - 760y^ + 1636?/^ - U17y^ + ASOy^ - 4y - 8))x^ 
+ (-16(9 - 6A)y^ + 8(1536^ - 1472)^^ + (43616 - 35586^)7/^ + (4251 6^ 
-109040)^/^ - 24(756^ - 5836)y^ - 8(36^ + 8428)t/2 + 96(6^ + 100) y + 192)a;^ 
+y(-2AbW + 24(136^ - 64)/ + (4672 - 12846^)/ + 33 (65^^ - 336)/ 
+ (42904 - 12906^)/ - 24(6^ + 2084)/ + 48(3 6^ + 268)y + 576)a;=^ 
+/(2462/ + (672 - 1866^)/ + 3(1616^ + 816) / - 4(1116^ + 1796)/ + 12(6^ 
-196)y^ + 32(36^ + UO)y + 576)x + y^ (-6b^y^ + (336^ - 80)/ - 6(96^ 
+ 196)/ + 12(6^ + 204)/ + 24 {b'^ - 28)y + 192)) 

+4(110/ - 293/ + 181y - 27) + 2(96/ - 522/ + 619y^ - 162y - 6)^^ 
-2y(6/ + 81y=^ - 282 / + U6y + 12)x + y^{6y^ + 33^^ - 76y - 12)) 

b^(-2y + xUy - 2) - 1) 

B^^iO) = ^ — ^ ^ ^ — / ' (B9) 

512v27r(x + y — 2xy)^Jx + y — 2xy 

for 7Jg'(x,y,bi,b2) and 

^10 (2) = 64(^^^-\)3(^ + ^)4 (-^'^' + + 9% - 64)x^ + (6^(9/ - 7y + 1) - 8(24/ 
+36y - 19)) x^ + ((16y^ - 267/' + lOy - 1)6' + 24(-24/ + 28^' + 14y - 3))x^ 
+(3(3 6' _ 64)/ + (800 - 266')y^ + 2(96^ - 200)^' - 3(6' + 8)y - 24) x' 
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+(96y^ - (76^ + 96)y^ + 2(56^ - 24)y^ - 3(6^ - 56)y^ - 64y + 24) x + y{-b^y^ 
+bY + - 40)i/-^ - (6' - 120)y' - 88y + 24)) 



+7)6^ + 24(8?/^ - 8?/ + 7)) + ((-16y^ + 58y^ - 50y + 1)6^ + 8(72y^ - 228y^ 
+222?/ - 79)) + {{-9y^ + 58y^ - lUy^ + 75y - 8)6^ + 8{2Ay^ - 212y^ 
+434y2 _ 301 y + 51))x2 + ((lly^ - 50y^ + 75y^ - 40y + 4)6^ + 8(-12y^ 
+142y=^ - 229 + 64y + 41))x + fe^y^ + (96 - 56^)y^ + (76^ - 152)y^ + (6^ 
-56)7/^ - 8 (6^ + 5)/ + 4(6^ + lU)y - 320) 



+20|/ - 3)x^ + 2|/(4 - 7y^ + 10|/ - l)a; - y^{y^ - 6y + 3)) 



A^f(2) = =((8y-7)x^ + (-24y^ + 30y-4)x^ 

32(0; + y - 1)^(2 - x - y)Wil ' - y) > K y y > 

+ {-2Ay^ + 122|/2 _ I46y + 39)3;^ + 2(4 + Iby^ - 73y^ + 71y - 13)x - 7y^ 
-Ay^ + 39y2 - 26y - 4) 

^ b'^{x'^ -Ayx + x + y"^ + y) ^ b\x^ - 4yx + x + y^ + y) 

' 64{x + y-l){x + yy ' ' 64{x + y - 1){2 - x - y)^' 



^r(2) = 7 777 T77 -(10b'^(l-2y)x^ + ((-72y'^ + 162y 

1287r(a; + t/)4(2-x-|/)4(a; + y- 2,T?j)2^ V VV y y 

-55)6^ + 128(1 - 2yy)x^ + ((192 y'^ - 210y + 49)6^ + 128(1601/^ - 256y^ 
+134y - 23))x^ + ((1160/ - 2992 y^ + 2240y^ - 765y + 150)6^ + 128(12y^ 
-584y=^ + 732y^ - 293y + 36)) x^ + (6^(1992y^ - 8460?/^ + 11538?/^ - 6591^^ 
+ 1868|/ - 228) - 128(320 y^ - 7647/^ - 290y^ + 7A7y'^ - 2My + 30))a;^ 
+ (6^(1608y^ - 9804 y^ + 19998y^ - 16821y^ + 6250y^ - 996y + 24) 
-I28y{<68y^ - 1004 y^ + 2142^^ - 799^2 _ 124?/ + 62))x^ + (128(56^ + 96)?/^ 
-16(341 6^ + 1600)/ + 2(79936^ - 35200)/ + (196992 - 195416^)/ + (9960 6^ 
-80896)^3 - 8(2476^ - 1728)^^ + 32(36^ - 224)y + 1536)a;3 + y(24 (56^ - 64)y^ 
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-96(156^ + 128)y^ + 32(1916^ + 656)y^ + (33920 - 11103 + (83626^ 
-96768)?/^ - 8(2796^ - 5824)?/^ + 144(6^ - 128) y + 4608)a;2 + y'^{l2h'^y'^ 
-6(296^ - 256)|/^ + 6(1736=^ + 640) y^ - (28456^ + 128)i/^ + 4(8516^ - 5248)7/^ 
+(33024 - 13806^)^2 + gg (^2 _ ^^^-^y ^ 4608)x - y^{<oh'^y^ + (384 - 576^)y^ 
+(2396^ + 384) y"" + (1536 - 4706^)y^ + 4(896^ - 1088)y^ - 8(36^ - bl2)y 
-1536)) 

~ 256v^7r(a; + yf{2 -x- yf{x + y- 2xyfl^ ~ ^ - (2?/ - l)(3(16y 

-5)6^ + 16(32y2 - 32y + 13)) - (36^(727/^ - 144|/^ + 647/^ - 13^/ + 3) 
-8(256y^ - 64y^ - 384 + 374y - 91) )x^ + (8(768y^ - 2688y^ + 1716y^ 
+786?/^ - 940?/ + 179) - 3 6^(128?/^ - 536?/^ + 650?/^ - 311?/^ + 74y - 8))x^ 
+ (8(256?/^ - 2688 + 7616|/^ - 5922|/^ + 1475|/2 - 297i/ + 124) - 36^(72/ 
-536y^ + 1124 - 833y^ + 221^^ - 8y - A))x'^ - (1024y^ + (512 - 4326^)^^ 
+50(39 - AU)y^ + (71056 - 24996^)y^ + 300(36^ - 176)y=^ + 16(36^ + 1165) y^ 
-48(6^ + 124)1/ + 864)a;^ + |/(246V - 96(6^ - lQ)y^ - 64(3 6^ - ?>2)y^ + 3(3116^ 
-7760)y^ + (56984 - 6636^)^^ - 48(6^ + 751) y^ + 24(36^ + 520)y - 2592)x^ 
+y^(-2Ab^y^ + (786^ - 672)y^ + (39 6^ - 3024)y^ + (13472 - 2226^)y^ 
+24(6^ - 995)2/^ + 16(36^ + 692) y - 2592)x + y^{<6h^y^ + (80 - 156^)2/^ - 9(6^ 
-104)y^ + 24(6^ - 111) y^ + 12(^2 + 296)y - 864)) 

Sr(2) = -r- — -r((10-20?/)x^ + (8y2 + 42|/- 15)a;^ 

^ ^ ' l2%Tx{x + yf{2-x-yf{x + y-2xyY^^ v y y ) 

+ (280|/3 - 436y2 + 218y - ^l)x^ + (72 y^ - 564y^ + 466^/^ - my + Q)x^ 
+y(12y^ - 102y3 + 378^' - 149y + 12) x + y'(-6y=^ + 33y' - 83y + 6)) 

^r(2) = ^ '- ((47/-2).x^ + (-I67/ 

256v27r(x + y){2 - x - y){x + y - 2xy)^Jx + y - 2xy 

+10y - 3)x' + (4y3 + lOy' - l)x - y(2y^ + 3 y + 1)) (BIO) 
ioTTg^{x,yMM) 
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APPENDIX C: THE FUNCTION s{x,b,Q) IN THE SUDAKOV FACTOR 



In this appendix we present the exphcit expression of the exponent s{x, b, Q) appearing 
in the Sudakov factor. Defining the variables, 

xQ , , 1 



In 



V2A 



6 = In ■ 



QCD 



bA 



QCD 



(CI) 



the exponent s{x, b, Q) is presented up to next-to-leading-log approximation [43| 



s{x, b, Q) 



g In ( X 



AO) 

\n{2q) + 1 ln(26) + 1 



+ 



In 



1 



^-1 



^4(2) ^(1) 



In 



In (I 



^3 [In2(2g)-ln2(26)] 



ln(2g) + 1 ln(26) + 1 



21n(26) + 1 



q b 

A^q-b 
432/?6 ^3 



16/?t 



21n(2g)+3 21n(26)+3 



9W(26) + 61n(26) +2 



18 ln2(2g) + 30 ln(2g) + 19 18 ln^(26) + 30 ln(26) + 19 



where the coefficients (3i and A^^^ are 

33 - 2n/ 



/5: 



12 



P2 



62 



153 - 19n/ 



(C2) 



24 

3 9 3 27 ^ 3 2 ' 



3 7b 



(C3) 



with 7^; the Euler constant. 

The exponent s{x,b,Q) is obtained under the condition that xQ / \^ > 1/6, i.e. the 
longitudinal momentum should be larger than the transverse momentum. So s{x, 6, Q) is 
defined for q > b, and set to zero for q < b. As a similar treatment, the complete Sudakov 
factor is set to unity, if e^^ > 1, in the numerical analysis. This corresponds to a 
truncation at large /ct, which spoils the on-shell requirement for the light valence quarks. 
The quark lines with large fc-p should be absorbed into the hard scattering amplitude, instead 
of the wave functions. 
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